Polarization transitions in interacting ring ID arrays 
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Periodic nanostructures can display the dynamics of arrays of atoms while enabling the tuning 
of interactions in ways not normally possible in Nature. We examine one dimensional arrays of a 
"synthetic atom," a one dimensional ring with a nearest neighbor Coulomb interaction. We consider 
the classical limit first, finding that the singly charged rings possess antiferroelectric order at low 
temperatures when the charge is discrete, but that they do not order when the charge is treated as 
a continuous classical fiuid. In the quantum limit Monte Carlo simulation suggests that the system 
undergoes a quantum phase transition as the interaction strength is increased. This is supported 
by mapping the system to the ID transverse field Ising model. Finally we examine the effect of 
magnetic fields. We find that a magnetic field can alter the electrostatic phase transition producing 
a ferroelectric groundstate, solely through its effect of shifting the eigenenergies of the quantum 
problem. 

PACS numbers; 73.21-b 



I. INTRODUCTION 

Fabrication of micro- and nano-sized structures such 
as quantum dots, wires and rings has made it possible 
for physicists to examine new ideas in electronic devices. 
These small size devices act as artificial atoms with spec- 
tra and shell structures similar to those of real atomsQ. 
However it is possible to control the properties of these 
synthetic atoms in a way that is impossible with real ones. 
For example, a regular atomic orbital is three dimensional 
and one has limited control over the electronic wavefunc- 
tion. In contrast, by controlling the shape of a quandum 
dot we can distort the wavefunction, controlling its po- 
larizability and its interaction with adjacent dots. While 
the properties of periodic arrays of atoms are well under- 
stood in solid state physics, we have new playing field - 
periodic arrays of nanostructures ~ in which we have an 
unprecedented control of the "atomic states." 

The focus of this paper is on the periodic one dimen- 
stional (ID) arrays of nano-rings. We chose this system 
for two reasons. First, quantum rings displayinteresting 
phenomena (e.g. persistent currents [J, |^ la, Q- ) which 
are not found in dots. The basic difference between the 
ring geometry and a quantum dot is the excluded middle 
which confines the electron in a ring to a narrow spatially 
periodic channel. This compact, periodic geometry can 
allow dynamics not found in other systems.]^ Second, 
it has been become possible to create extremely small 
rings. These arrays of nano-rings can be fabricated ei- 
ther by dry etching P| or by using MBE techniques to 
foster self assembled InGaAs/GaAs rings. The size of 
these nano-rings is ^30 nm for outer radius and ^10 nm 
for inner radius for self-assembled InGaAs/GaAs rings 
0- Such techniques not only produce extremely small 
rings, they also make it easy to make periodic arrays of 
small rings. 

In this paper we consider an ideal array of ID rings 
at zero temperature, each carrying a single charge. 
The rings are sufficiently close together that there is a 



FIG. 1: A schematic picture of the groundstate of classical 
point electrons for ID array of rings. The ring radius is R 
and the separation is D. The ID ordering is antiferroelectric 
for the infinite size system and thus has a double degenerate 
groundstate. 



Coulomb interaction between the electrons, but sepa- 
rated enough so that the tunneling between rings can be 
neglected, [3 as discussed in section Ull In section Hill we 
consider the classical case in section and see how different 
assumptions allow for symmetric or symmetry breaking 
ground states. Section Hvl contains the main results of 
this work, where we show that there is a quantum phase 
transition in a ID array of rings for _B = 0. We show this 
both through Monte Carlo simulation of 1 -|- ID classi- 
cal statistical respresentation of the problem as well as 
mapping it on to the ID transverse field Ising model. 
The polarization pattern is antiferroelectric. In section 
0we examine how this transition is affected by magnetic 
fields. We conclude in section|^by summarizing results 
and discussing possible applications. 



II. THE MODEL 

We consider a one dimensional array of singly charged 
narrow quantum rings with radius R and center-to-center 
separation of D (FigQ. The width of each ring is much 
smaller than its inner radius so that we need only consider 
the one dimensional movement of the electron around 
the ring. While the rings are isolated from each other so 
there is no charge transfer between rings, the rings still 
interact electrostatically with their nearest neighbors. (9| 
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FIG. 2: Normal modes of the ID ring array with gaps of (solid 
line, inset (a)) 2^/2ujo and (dashed line, inset(b)) 2\/Qijjo, 
where ujq — /m*D'^ 



While in principle the Coulomb interaction is long range, 
we assume that there is sufficient screening that next 
nearest neighbor interactions can be neglected. 

We also neglect tunneling between rings. This is rea- 
sonable since the physical realizations discussed above 
produce rings that are not in tunneling contact. Tunnel- 
ing will be exponentially suppressed with D, so that for 
spatially distinct rings it will be minimal. 

All the phenomena explained in this article only will 
appear in experiments if the electrons do not lose their 
quantum mechanical phase, i.e. the ring's perimeter 
has to be smaller than the electron's coherence length 
{2ttR < ) and the temperature has to be lower than 
the dephasing temperature (T < T^). 

In each ring the confinement energy of the electron 
scales as Eq = l1m*T(^\ this energy opposes local- 
ization of the wavefunction in the ring. The inter-ring 
Coulomb repulsion, which scales as = /D, tries to 
localize the wavefunction. Electrons are repelled from re- 
gions of the ring where it is too close to the charges on 
neighboring rings. The competition between these two 
physical scales creates a quantum phase transition in the 
array from a localized state to extended state as we will 
see below. 



a continuous fluid. 



A. Classical point charges 

The classical model considers one charged point parti- 
cle per ring with only nearest neighbor Coulomb interac- 
tion. Unlike the quantum mechanical case there is only 
one energy scale in the classical problem which is the 
Coulomb energy Ec — /D. The energy of a ID array 

is given by Um = Ya^i /\n{Oi) - ri+i{Ot+i)\, where 9i 
is the location of the i-th electron as measured from hor- 
izontal axis. In the dipole approximation we can write 
this as : 
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Si+i + - {D ■ (si- Si+i) f 



(1) 



where e = R/D and Uq is a constant, Uq = + y)- 

In the second expression we identify the position of each 
charge by a vector Sj in the 2D plane pointing from the 
center of the i-th ring to the charge on that ring. The 
unit vector D lies on the horizontal axis. 

The cos 20 (or {D ■ s)^) term explicitly breaks the ro- 
tational symmetry, driving the system from XY to Ising- 
like behavior. The Heisenberg term in the last line of 
equation drives the system ferroelectric at zero tem- 
perature while the second and larger term favors states 
where neighbors point in opposite direction. Thus the 
system at zero temperature orders in an antiferroelectric 
(AFE) pattern (FigQ in one dimension. Our numerical 
Monte Carlo simulations of the exact Coulomb interac- 
tion also verifies the existence of such a minimum energy 
configuration in the classical finite size arrays. 

We can examine the stability of the AFE state by 
finding the higher energy modes of the system. We 
expand the energy function ^ to quadratic order in 
displacement angle around the AFE configuration using 
di = {—ly^+ai. The AFE configuration has a basis 
with two sites so we find two independent normal modes 
with frequencies : 



III. CLASSICAL RESULTS 
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Before solving a quantum mechanical problem it is of- 
ten helpful to look at the similar classical case which is 
usually easier to solve. Below we consider two classical 
models, one in which the classical charges are treated as 
ideal points, and the second in which they are treated as 



Where ujq = y/e^Jm*D^. Both the modes are gapped 
since the Ising-like term provides the harmonic restor- 
ing force at each site. The modes are shown in Fig.||2Jl. 
Normal modes are found to be independent of the ring 
radius. 



3 




0.1 0.2 0.3 0.4 O.E 
R/D 



FIG. 3: A plot of the second fourier amplitude of the classical 
charge distribution on a ring in a ID horizontal array. The 
circles are numerical results, the solid line is a scaled plot of 
eq.Q. Scaling is required since the analytic result neglects 
all higher fourier modes. Inset: a sketch of the charge distri- 
bution that corresponds to this Fourier mode. Note that the 
symmetry of the array is not broken by the charge distribu- 
tion. 



minimum energy configuration of the ID array of charge 
fluid does not break the up-down symmetry of the sys- 
tem. 



IV. QUANTUM RESULTS FOR B ^ 0. 

At first glance the quantum mechanical wavefunc- 
tion of a charged particle resembles the classical charge 
fluid. Although the wavefunction does not have a self- 
interaction, the quantum particle has a kinetic energy 
which opposes localization, making the analogy to charge 
fluid even stronger. Unlike the classical case, the quan- 
tum problem has two competing energy scales: the 
quantum kinetic energy. Eg, preventing localization and 
the Coulomb interaction energy, Ec, trying to force the 
charges away from each other. At Eg ^ Ec we expect 
charge localization on each ring and at Eq 3> Ec we ex- 
pect no localization of charge. However it is not a priori 
obvious whether the charge localization (system polariza- 
tion) breaks symmetry or not, whether this localization 
is a smooth function of the external parameters, and if 
it is a phase transition, what is the exact nature of this 
transition. 



B. Classical charge fluid 

Another interesting classical limit of our ring problem 
is when there is a classical self-interacting fluid of charge 
on each ring while the nearest neighbor fluids are still 
interacting with each other. 

To find the minmum energy distribution of charge den- 
sity on each ring we define an angular dependent charge 
density Pi{Oi) on each ring where: J pi{9i)d9i = 1. We 
are looking for the minimal solution to the variational 
quantity : 
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de 



XJ2 d9pM- (3) 



For a ID ring this expression is divergent due to self 
energy. We can regularize this in several ways. One 
method is to introduce a short distance cutoff C to the 
Coulomb interaction, discretize the integral equation and 
then solve the problem numerically. An approximate an- 
alytic solution can then be obtained by Fourier expanding 
the distribution, keeping only the first three modes. 

For ID array with periodic boundary condition we find 
that the amplitude of the non-trivial Fourier mode as a 
f miction of ( and e = R/ D ior p is given by : 



n(2) 



-3^£3(2-5C') 
4(-2 + 41og(e/C))' 



(4) 



We compare this analytic result with the numerical diag- 
onalization of Q in Fig.©. As we can see in Fig.© the 



A. Variational calculation 

As a first step we can use a simple variational wave- 
function to find the polarization behavior of the system 
in ground state. The dimensionless hamiltonian of an 
array of one dimensional rings with radius R is given by: 
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\rM)-rA0i)\' 



(5) 



where S — Ec/Eq is the interaction strength and the 
energy is measured in units of Eq = /2mR^. 

To find the ground state energy of the ID array we 
employ a simple ansatz for the wavefunction of each sub- 



lattice : iPa{&) 



y- 



/27r 



3(61 - (f)) and t/jsiO) 



cos{0 — (j)) alternately. The ground state 
is obtained by minimizing the energy 
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values of y and 

using dipole approximation for Coulomb interaction 
we find: 
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for 8 > (5c (e) 
S < Sc{e) 



(6) 



and (j) = it/2, where the critical value of the interaction 
is given by Sc{e) = To find out the degree of 

polarization we define the staggered polarization vector 
as: 



(7) 
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FIG. 4: A comparison of numerical and analytical calcula- 
tions of the staggered polarization and energy as a function 
of 5 in a ID quantum ring array obtained in the Hartree ap- 
proximation. The numerical results are for the case include 
Fourier modes |m| < 1 (triangles) and |m| < 6 (boxes). The 
solid line is the analytic result assum- ing \m\ < 1. The quan- 
tity (5 is a measure of the competition between the Coulomb 
interaction and the quantum kinetic energy 

Using variational results the staggered polarization of the 
system is as follows: 

1 S<Sc 

Where D± is the unit vector perpendicular to the com- 
mon axis of the rings (Fig^ and 5c is defined in equation 

As we can see variational calculation suggests that the 
ground state of the ID array of rings antiferroelecrically 
polarizes in perpendicular direction at high interaction 
strengthes while at lower values the wavefunctions are 
not localized, hence the system has no polarization. The 
validity of this result will be confirmed in next sections 
using more exact and reliable methods of calculation. 

B. Hartree Approximation 

The rings considered here are well separated with ex- 
actly one electron on each ring. Under this condition and 
because of strong Coulomb repulsion the effect of inter- 
ring transfer of electrons and overlap of wavefunctions is 
small. We can therefore neglect the inter-ring transfer 
from our calculations. Since without overlap the elec- 
trons do not have any exchange interaction, the Hartree 
approximation is exact for this problem. [Tol| 

We can decompose the wavefunction in each ring 
into a limited number of Fourier modes, ipii^^) = 
Sn°=-n CnC™*, and then solve the system numerically 



in the Hartree approximation. We impose the periodic 
boundary conditions on the array and by an iterative self- 
consistent method we find the ground state wavefunction 
of the rings. In Fig.Q) we can see the numerical results 
of the polarization and energy change of the ID array 
of rings for different number of Fourier modes using ex- 
act Coulomb interaction and also its agreement with the 
variational calculation when we restrict the number of 
Fourier modes to n E { — 1,0,1}. The results are little 
changed when we increase the number of Fourier modes 
mostly in the high coupling regime. 

All the above results suggest that there is transition 
from unpolarized to polarized state at zero temperature 
by changing the coupling. By looking at the behavior of 
polarization when the number of Fourier modes increases 
we realize that this transition tends to be sharper and 
sharper for higher number of Fourier modes suggesting 
a true phase transition in the system. If true, this tran- 
sition would be a sudden change of ground state of the 
quantum system at zero temperature, known as quantum 
phase transition. [Tlj We demonstrate that this is the case 
and determine the universality class of the transition in 
the next section using the Monte Carlo simulation. 



C. Monte Carlo simulation 

It is well-known that we can write the quantum par- 
tition function of a quantum system, Z = Tre^^^ as 
the sum over all paths taken by the system in imaginary 
time defined by the scale ?l/3. If the quantum system 
is D-dimensional then the partition function will look 
like the path integral of a D-|-l-dimensional classical sys- 
tem in which the extra dimension is the time direction 
< r < pti. At zero temperature (3 ^ oo the clas- 
sical system is truly D-|-l-dimensional. One can derive 
an effective Hamiltonian for such a classical system from 
the quantum Hamiltonian using a complete set of ba- 
sis states. In this classical system the parameters of the 
quantum system (in our case 5) is a control knob like tem- 
perature. We can use Monte Carlo simulation of such a 
classical system and find out the universal behavior of 
the quantum system. 

To develop a l-|-l-dimensional classical theory for our 
ID ring array we first stagger the order parameter, 9i — > 
{—lyOi so that we can analyze the Monte Carlo results 
easily. We also use dipole approximation for the Coulomb 
interaction. Consequently we can write the Hamiltonian 
of the system as: 

j=i ^ j=i 

In which = h^mR^ and Ej = e^t^j^D. The stan- 
dard derivation|l2| using the Villain approximation [l^ 
IT^ tells us that the 1-l-lD classical partition function 
equivalent to the ID interacting quantum ring array at 



zero temperature is: 

AT 
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fc=l 



OkiTa)] 



(10) 



Where V{0) = De{Ta) and 

yfc(Ta) = 3cOs[0k{Ta) - 9k+l{Ta)] 



7fe+l 



C0s[6'fc(Ta 

- 3cOs[9k{Ta) - 9k+l{Ta)] C0s[9k{Ta) + Ok+l{Ta)]. 

(11) 

The parameter r has the dimension of time and AfSr = 
f3h. It can be shown that the field 9{x,t) obeys the 
periodic boundary condition , 9{t + (31%) = 9{t) (^3)- 
We will also assume periodic boundary condition in space 
direction all over the simulation. 

By defining the spin vector Si — (cos0i,sin0i) we 
can interpret Eq. H1U|) as a two dimensional classical spin 
model. Our early calculations suggested that the sys- 
tem of ID rings has a transition from the unpolarized 
to the AFE state. In this classical analogue because we 
have already staggered the order parameter we expect 
to see a transition from unpolarized to ferromagnetically 
polarized state (FE). Close to this transition the spatial 
variation of the order parameter 5* is smooth so we can 
approximate (|ll|l as follows: 

Vk{Ta) KiC0s{9k{Ta) - 0fc+l(T,)) - 2 COS 2^^fe (t,,) . (12) 

Using the above potential finally the classical partition 
function looks like: 

Z (x D9{t) exp {k^cos(0. -(?,)- —^cos2(?,}. 

(13) 

Where i and j run over an infinite 2D square lattice and 
we have determined 6t to identify the two couplings in 

ea. (P|l as if = • Equation (O is a 2D XY model 

with a symmetry breaking field which is 2/3 the XY cou- 
pling. Our Monte Carlo analysis shows that this model 
has a continuous phase transition. The order parameter 
of this system is the total magnetization density equiv- 
alent to the total staggered polarization of the ID ring 
array: 



P. 



(14) 



Where the average on the left hand side is the thermo- 
dynamic average over the infinite size lattice. 

The fluctuations in this system is controlled by K 
which is the analogue of 1/T in real classical systems. 
We can measure the analogue of specific heat of the sys- 
tem using: 




FIG. 5: Monte Carlo results of the Cv = AE'^/N'^ for different 
system sizes. The system is a 1-l-lD classical equivalent of ID 
quantum ring array at zero temperature. 
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FIG. 6: Plot of critical coupling Kc{L) at different system 
sizes taken from the Cv plots. The solid line is a linear fit to 
the data indicating Kc{oo) « 0.699. 



in which (.) is the average over an ensemble and E is the 
total energy of the Nx N system. This quantity diverges 
at the critical point of the infinite system undergoing a 
continuous phase transition. Fig.© shows the change of 
the specific heat of our 1-1- ID system in terms of the pa- 
rameter K for different lattice sizes. As we can see at Kc 
the peak gets sharper and sharper with increasing lattice 
size L. An extrapolation of the point of the maximum 
of Cv, Kc{L) to = determines the approximate 
critical point of the infinite lattice(Fig|Bl) . Also an ex- 
trapolation oim{L, K) for different values of K in Fig.© 
shows that a real continues phase transition happens in 
the infinite size system. 

The effective classical system derived here does not 
fully explain all the physical aspects of the ID quan- 
tum system mainly because of the approximations used 
to derive the path integral. However, we believe that, 
close to the critical region, these approximations do not 




of the 1+lD classical system to infinite size at couplings -f^ 
i<'/ii'c(oo) = 1.45(empty boxes), 1.32(triangles), 1.03(poly- 
gons), 0.74(filled boxes), 0.45(stars) and 0.16(crosses). The 

solid fines are linear fit to each set of data. FIG. 8: Plots of Binder ratio for diflterent system sizes. The 

behavior is sharper at larger sizes. 



play any role in the general behavior of the system and 
the universality class remains unchanged. Hence using 
the finite size scaling method we can determine the crit- 
ical exponents of the classical system and determine the 
universality class of the actual quantum system. 

1. Finite size scaling of the 1 + lD system 

One of the best parameters for examining the phase 
transition and find the universal exponents with finite 
size scaling is the dimensionless Binder ratio: 



defined for a system with size L. In the disordered phase 
K < Kc the correlation length ^ is finite so for L ^ ^ 
the distribution of m is Gaussian around m = with the 
width A^"i/2 - L-'^/^ so g ^ 0. On the other hand 
for K > Kc where (m) is finite, approaches a constant 
as L — > oo. The variation of with K becomes sharper 
and sharper as L increases, however all the g's cross at 
the transition point Kc- The variation is given by the 
following finite size scaling function: 

gL{K)^g{L^/^{K-Kc)). (17) 

Where ^ is a scaling function which depends on L and 
K only in that particular form. By using the finite size 
data we can try to find a data collapse and by calculating 
the standard deviation find the best exponent v fitting to 
the collapsed function. Fig® shows the Binder ratio for 
different lattice sizes. Fig.® shows the collapsed data 
and Fig.l^UJ shows the best exponent is v — 0.99 ± 0.01. 
The error is estimated from the mesh of the numerical 
calculation. 

The scaling for the order parameter |m| is: 

m = L-f^/^'X^iL^/'^iK ~ Kc)). (18) 



3 






2.5 




■ N=4 






N=8 


2 




- - N=10 






— N=20 


1.5 
1 







-0.5 0.5 1 

FIG. 9: Collapse of Binder plots at the critical region. The 
best collapse is obtained for p — 1.01 ± .01. 

Where is a function of a; = L^^'^{K — Kc) only and 
(3 is one of the universal scaling exponents of the system. 
To determine the universal exponent (3 we plot L^/" vs. 
X for different sizes. Fig. H12|l and H13I) show the the col- 
lapse of different data sets and the standard deviation 
for different exponents respectively which shows the best 
estimation is /? = 0.125 ± .005, or 1/8. 



2. Universality Class 

The universal exponents extracted from the finite size 
data indicate that our 1-1- ID classical XY model in the 
symmetry breaking field is in the universality class of 2D 
classical Ising model, hence the nature of quantum phase 
transition of our ID quantum ring array is Ising-like. The 
Coulomb repulsion forces the electrons to alternate stay- 
ing on the top and bottom of the rings. However the 
quantum kinetic energy tries to avoid localization. This 
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FIG. 10: Standard deviation of the set of scaled plots of 
Binder ratio for different exponents. The case for v = 1 is 
the best choice which is plotted in Fig.©. 



FIG. 12: Collapse of scaled magnetization data sets for v — 1 
and /3 = 1/8 in the critical region. 
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FIG. 11: Monte Carlo results of average magnetization of the 
1+lD classical system for different system sizes. The change 
of magnetization tends to be sharper as the system size grows. 



FIG. 13: Standard deviation of the set of scaled plots of mag- 
netization for fixed v = 1 and different exponent 13. The case 
for /3 = 1/8 is the best choice which is plotted in Fig. 1121 . 



kinetic energy causes the electrons to tunnel from top into 
bottom of the ring hence destroying the antiferroelectric 
order. This ordering behavior shows up in the probability 
distribution on each ring. Fig. (|14|l shows the energy of 
each electron with the wavefimction ipdiO) = x+y cos(20) 
compared to when the wavefunction is a constant all 
around the ring. The wavefunction has two maxima 
on top and bottom of the ring which means the electron 
is fluctuating up and down. As we can see by increasing 
the coupling the lower energy state selected by the exact 
Hartree calculations (dots) gradually matches ipd instead 
of the constant wavefunction. This behavior persists in 
a range of couplings close but smaller than the critical 
couphng i.e. in disordered region 5 < 5c- Needless to say 
that after transition point the ground state wavefunction 
is no longer ipd and the system starts to excite more an- 
gular momentum eigenstates (Fourier modes). 

All the above discussion suggests that the nature of the 
antiferroelectric transition is not just the simple 2D Ising 
but is similar to ID transverse field Ising{T¥l) which 
has a quantum phase transition at zero temperature in 
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FIG. 14: Plot of the energy of interacting quantum ring ar- 
ray when all the wavefunctions are constant around the ring 
(dashed) or all are in the form oi ipd = a + 6cos 20(solid). 
The points are the actual results coming out of the numerical 
Hartree calculation indicating that is the selected behavior 
for 5 < 5c ~ 8. 
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FIG. 15: Plot of the cos 26 potential around a ring (the thick 
solid line) , the ground and first excited state of this potential 
(thin solid lines) coming out of a simple numerical Schrodinger 
equation solver and the up and down states (dashed lines) con- 
structed from the two eigenstates(see appendix). The scale of 
the potential is exaggerated for easier comparison. 



the same universality class as 2D Ising. We can develop 
an effective ID TFI hamiltonian for our ring array in 
the dipole approximation. In this approximation we can 
write down the hamiltonian(|51) as follows : 



H = Ha + V 
Ho 



E 



HTTTrf + iSe^ cos 20, 

Cltli 



V = -,5e2^[3cos(0, + 0,+i)+cos(0,-0,+i)] 



(19) 



The cos 29 term in hamiltonian Hq has two minima at 
top and bottom of the ring. Fig. (|15|l shows the poten- 
tial and the two lowest energy states of it with energies 
Eq < El. The rate of tunnelling from top to bottom or 
vice versa is determined hy A = Ei — Eq. The poten- 
tial V tries to align the electrons hence it acts like the 
Ising interaction. A more rig orous derivation using the 
Holstein-Primakov bosons shows that the projection 
of the hamiltonian Ti into the subspace of the ground 
and first excited states of Hq can be written as (see ap- 
pendix) : 



N 



H 



N 



'i+i- 



(20) 



iSe''. Nu- 
l-0.1(5e2 



In which cr's are Pauli spin matrices and J 
merical diagonalization of Hq tells us that A 
for small e. 

Close to transition the Coulomb repulsion is not strong 
enough to excite the electrons to higher states, conse- 
quently the TFI model in Ea. (|20|) is valid and indicated 
the nature of transition of the ID ring array. 



V. QUANTUM RINGS TRANSITIONS FOR 

B / 

Y. Aharanov and D. Bohm (AB) have predicted that 
the wavefunction of an electron moving in a vector po- 
tential A{x) along the path C acquires a phase shift: 



AA 



e 

he 



A ■ dr. 



(21) 



AB predicted that this phase shift can be observable. 
When an electron is confined on a closed path like the 
case of charged ring threaded by magnetic flux (f> the 
phase shift after one 27r rotation would be : 



he 



A ■ dr 



(22) 



Where (j>Q = he/e ~ 4.135 x 10~^G.em^ is the quantum 
of flux. The phase shift above has been observed in nu- 
merous experiments and different devices including the 
experiments of persistent current and excitons in quan- 
tum rings. In this section we show how magnetic field 
changes the behavior of polarization. 

The Hamiltonian of an electron in ID a ring threaded 
by a constant uniform magnetic field Bz (the ring is in 
the x-y plane) is: 



(23) 



in which the choice of gauge: A — -|-(— ?/, x, 0), the mo- 



mentum is in polar coordinates:^ 
eigenfunctions are periodic: ^p{9 + 
eigenenergies of will be: 



27r) 



'Roe 

= m- 



En = 



2TOi?2 



The 



(24) 



in which n is an integer. By changing the gauge A A — 
VA wavefunctions undergo a phase change tp ~* e'^^ip. 
For example we can use the gauge transformation with 
the choice of A = [BE? / 2)9 to remove the vector poten- 
tial from (|23|l but at the same time we have to shift the 

phase of the wavefunctions to e~^^^^. As the result the 
eigenfunctions change: 



M0) 



1 



J(n-4>/(po)6 



2tt 



(25) 



This eigenfunction however has a different boundary con- 
dition than the previous one : 



^{9 + 2tt) = e2"'^/'^''V(6') 



(26) 



but as one physically expects the eigenenergies are not 
changed. The ground state energy of the interacting 
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FIG. 16: Results of numerical Hartree calculations of the po- 
larization of a ID quantum ring array threaded by half-flux 
quantum (p/cpo = 1/2. For S < Sc the system displays a lon- 
gitudinal ferroelectric (FE) polarization (squares), while for 
5 > 5c an antiferroelectric (AFE) polarization (triangles) is 
observed. 
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FIG. 17: Plots of ground state energy of the interacting quan- 
tum ring array in the external magnetic flux threading each 
ring for different couplings 5. The physics is periodic because 
of Aharanov-Bohm induced phase that is proportional to the 
flux. 



quantum ring array in the magnetic flux can be written 
as: 



n— — no 



\m-r,i0') 



(27) 



where *o(^) 



J2n=-no Cne" IS the ground 



state wavefunction expanded in the free hamiltonian ba- 
sis states. As we can see in Ea. (|27|l the only part that 
is affected by the AB phase is the kinetic energy and the 
potential energy is not sensitive to the phase. The en- 
ergy calculated in Ea. (|27(l is periodic in (po since when 
(f) = (jjQ = I we can rearrange the infinite sum and show 
that the value of the kinetic energy is equal to its value at 




2 7T 



2 TT 



FIG. 18: Plot of the two-body potential function in Eg 1281 
The interacting rings in magnetic field select the minimum 
of this potential for their ground state wavefunction at low 
couplings. The potential is in units of 5e^. 



(j) — 0. This is due to the well-known fact that the physics 
of quantum rings does not change at integer flux quanta. 
For infinite uq this argument is true at any range of mag- 
netic flux; in our numerical calculations where we have 
used a finite number of Fourier modes Eq is periodic only 
in a finite range approximately given by < < no0o in 
Fig. (|17|l we can see the periodic behavior of the ground 
state energy of the ring array as magnetic flux changes. 

The results of our numerical Hartree calculations in- 
dicate that in a ID ring array in which each ring is 
threaded by a magnetic flux cj) the polarization pat- 
tern changes from unpolarized to ferroelectric at half- 
integer flux quantum. Fig. (|16|l shows the behavior of 
Px , the component of the total polarization vector 
P = J2, de\'tP,{e)\^fi{e) in the direction of the ring's 
common axis at half flux quantum. This plot shows that 
at (5 < (5c the wavefunction has an unbalanced distribu- 
tion around each ring. However the total polarization 
vanishes at higher values of interaction where the wave- 
function distribution becomes antiferroelectrially polar- 
ized in the array. The finite polarization at small in- 
teraction strengths has a ferroelectric pattern which is 
degenerate left or right. From this result we can see that 
the physics of quantum ring arrays changes at half in- 
teger flux quantum. The total staggered polarization in 
the y direction perpendicular to the common axis of the 
rings starts to build up at (5 > dc as in the case of no 
magnetic field. We can explain this phenomena of finite 
transverse polarization due to magnetic field in different 
approaches. 

We can use a simple perturbative discussion to under- 
stand this behavior qualitatively. The Ea. (|19|) which is 
the dipole approximation of the total hamiltonian will 
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modify in presence of a magnetic flux as follows: 



H = Hn + V 



Ho 
V 



E 



-36e^ cos 29,] 



-5e^ J2 [3 cos(0, + + cos(0, - 



(28) 



In the above equation the kinetic energy hamiltonian has 
a degenerate ground state. For example at half filling, 
0/(/)o = ^,n = O and n — 1 levels are degenerate un- 
like the case of zero magnetic field in which the ground 
state is unique and at n = 0. By adding the symmetry 
breaking term cos 26 in the zero flux the electron gains 
enough energy to excite to the next higher level. This 
excitation causes the electron to destroy any localization 
in the range where V is not strong enough yet. However 
when there is a finite magnetic field the cos 29 can not 
lift the degeneracy between n = and n = 1 levels in the 
range of small (5's. That is why in this case the ground 
state of Hq remains degeneratefFig lTTIl . As long as 6 is 
small the perturbative two body potential in (|28|l can not 
excite the electron to higher levels and the kinetic energy 
of the electron freezes. When this happens the electrons 
behave classicaly and choose a wavefunction that mini- 
mizes the potential In Fig. lfTF|) we can see a 3D plot of 
the two body potential V in which it has two stable min- 
ima at (tt, tt) and (0, 0) indicating the preferred state of 
the quantum ring array at low S being the ferreoelectric 
right or left state. 

It is surprising that the system orders ferroelectricalUy 
along the chains. This occurs because the charge distri- 
bution is quite broad. Using only the two lowest Fourier 
modes the charge distribution can be written in gen- 
eral as p{9) = (1 + cos {9 — 9o)) /27r , which nearly wraps 
around the ring. If we change this to an (unphysical) flat 
charge distribution with an artiflcially varying width, the 
FE state has a lower energy than the AFE state when the 
width A9 ^ Abn, while for a triangular distribution the 
FE distribution is favored when A9 ^ O.Gtt. While the 
exact transition depends of course upon e, these calcu- 
lations help explain why the FE phase wins out at half 
filling. 



VI. CONCLUSIONS 

In this paper we have shown that there is a phase 
transition in a periodic array of electrons each confined 
to a ID ring. The parameter d determines if the ar- 
ray will spontaneously polarize; in ID the transition is 
at (5 « 10. It is easy to achieve small values of 6 sim- 
ply by choosing the ring separation to be large. Thus 
the "quantum" limit where the kinetic energy dominates 
is simple to obtain. To obtain the antiferroelectrically 
ordered state we need large S. We may write this as 



6 = {B? /aoD) X {m* /m) where qq is the Bohr radius 
and TO* is the effective mass of the electron. If R and D 
are R and D measured in nanometers, and to = (m* /m), 
then S ~ 18.9(i?^/£>)TO. We require that the rings do not 
intersect, so that D > 2R. Thus the ability to achieve 
large values of S in semiconductors will depend upon the 
value of the effective mass. If we set D = 2R, then for 
GaAs (m — 0.06) ID arrays of rings with a radius greater 
than ~10nm will be polarized. For AlAs (m — 0.4) the 
crossover radius is about 70nm. Rings with a smaller 
radius will not spontaneously polarize, but instead be 
isotropic 

It is well known that in ID there is no ordered state 
for T > for the Ising model. However, for small ar- 
rays over finite time intervals the system can order. To 
observe this behavior we want the characteristic energies 
of the system to be greater than the temperature. For 
the Coulomb energy kT < /D, which we may write 
as i)T < 1.8 X 10^ where T is in Kelvin. For the ki- 
netic energy this means kT < fi/2m*R^; if we measure 
{m/m*)R'^T > 40 in the same units. For GaAs we can 
choose R to be about 14nm at 4K; choosing materials 
with a smaller effective mass or going to lower tempera- 
ture allows us to increase the radius. 

An AFE polarized ring array will scatter light at a 
wavelength commensurate with the inter-ring separation, 
D. In ID there is a gap v^wo, which we may write 
as 2\/2rh{ao/ D)^^^ . For GaAs rings with a separation 
D = lOOOnm this gives w - 6.0 x lO^^Rz. The 2D arrays 
have a similar sized gap at zone center, but the gap van- 
ishes at one zone edge. The excitation spectrum can be 
probed optically, but scattering at the edge of the zone is 
difficult due to the constraints imposed by conservation 
of energy and momentum. Typically in such cases Ra- 
man scattering can be used to investigate the excitations. 

While we have not explicitly addressed the 2D case 
here, much can be gleaned from our results. The 2D clas- 
sical problem obviously has a finite temperature phase 
transition, as shown by our Monte Carlo simulations. 
The 2D quantum problem can be mapped on to the 3D 
XY model, which is known to order. We have performed 
simulations on the 2D case and find that it orders in a 
striped phase. 

Finally these calculations assume that each ring is 
singly occupied. This might be obtained by fabricating 
the rings upon a thin insulating layer covering a gate. By 
tuning the gate voltage we can bias the system so that it 
is energetically favorable for an electron to tunnel to the 
rings. The gate will also serve to cutoff long distance in- 
teractions between the rings, supporting the assumption 
of the nearest neighbor interactions used here. Moreover, 
this letter serves to start investigation into a broad class 
of problems, such as rings occupied by an optically ex- 
cited exciton/hole pair or perhaps by a small, varying 
number of electrons created by a random distribution of 
dopants. 

The topic of quantum dot arrays and their correlations 
has obvious and useful analogies with solid state models 
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of crystalline arrays of atoms. In this paper we wish to 
point out that experimentalists have at their disposal a 
host of "unnatural atoms" analogs: rings, quantum dot 
quantum wells, quantum rice, etc. The electrons in these 
nanoscale constituents are confined to orbitals that may 
not have atomic analogs. Morever, it may be possible to 
tune the shape of the constiutent to optimize some de- 
sired collective property such as frustration in electric or 
magnetic polarization, high susceptibility or sensitivity 
to optical polarization of light. Even more rich behavior 
will develop if we allow electrons to tunnel between these 
nanoscale periodic structures. 
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APPENDIX: TRANSVERSE FIELD ISING 
MODEL 

In this section we explain how one can write the projec- 
tion of the hamiltonian l|19|) into ground and first excited 
state subspace of Hq as a ID transverse field Ising model 
fEa l20|l . In order to make the analysis easier we change 
the variables 9i to staggered one, {—lyOi. Defining |0) 
and |1) as ground and first excited states of Hq and Eq 
and El the corresponding eigenenergies we can write the 
up and down states fFig ll5|) as: 

1 



So we can write the hamiltonian (|19|l as: 



IT) = ^(|0) + |i)) 
li) = ^(|i)-|o)). 



(A.l) 



Now we define the creation and annihilation operators: 

4it> = 4u) = o 

4U) = IT) , 4lT>-U), (A.2) 

and one can show that : 



N 



i=l Q,/3=T,i {ij) klmn=T,l 

(A.4) 

in which i,j are spatial indices and: 



1 so A 

2 I A £0 



(A.5) 



where eq — Ei + Eq and A — Ei — Eo- To calculate ma- 
trix elements of the potential we use a simple numerical 
Schrodinger equation solver to find the following quanti- 
ties: 

(T |cos0| T) = (i \cose\ i) «o 
(T |cos0| i) = (i \cose\ T) «o 

(T Isin^l t)«+l , (i|sin(?U)« 

(T|sin0U) = a|sin0|T)«O (A.6) 

In the above diagonal matrix elements of cos^ are ap- 
proximately zero because must of the wavefunction is lo- 
calized around 9 = ±7r/2. For the same reason the diag- 
onal matrix elements of sin 9 are ±1. Using the above we 
can derive the potential in the following second quantized 
form: 



V = -2<5e2^(4^c]^a 

iv) 

~ ^li 4t S"T + ^li 4i S"i ) 



t J 



(A.7) 



(A.3) 



Now we can use the Holstein-Primakov transformation 
[l^ to construct the following SU(2) covariant operators 
for each lattice point: 

5*+ = c\c^ , S- = c\c^ , S'^ = i(cjc^ - c\c{). (A.8) 

Using the above definitions and the fact that: cjc| -|- 
cjcj 1 for each lattice point we arrive at the following 
expression for the hamiltonian H19(l : 

N N 

J? = -8<5e2^5f5f+i+A^5r, (A.9) 

i=l i=l 

in which: S^'^y = {S+ ± iS-)/2. 
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